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A bstract 

This paper describes the interaction of some of the important 
parameters affecting the mixing process in a quick mixing region of a 
rich bum/quick mix/lean bum (RQL) combustor. The performance 
of the quick mixing region is significantly affected by the geometric 
designs of both the mixing domain and the jet inlet orifices. In this 
paper, several of the important geometric parameters and operating 
conditions affecting the mixing process were analytically studied. 
Parameters such as jet-to-mainstream momentum flux ratio (J), mass 
flow ratio (MR), orifice geometry, orifice orientation, and number of 
orifices/row (equally spaced) around the circumferential direction 
have been analyzed. Three different sets of orifice shapes were 
investigated: (1) square, (2) elongated slots, and (3) equilateral 
triangles. Based on the analytical results, the best mixing 
configuration depends significantly on the penetration depth of the 
jet to prevent the hot mainstream flow from being entrained behind 
the orifice. The structure in a circular mixing section is highly 
weighted toward the outer wall and any mixing structure affecting 
this area significantly affects the overall results. The increase in the 
number of orifices per row increases the mixing performance. 

Strong swirl in the quick mixing region appears to have an adverse 
effect on mixing at higher J conditions. Higher slot slant angles and 
aspect ratios are generally the best mixing configurations at higher 
momentum-flux ratio (J) conditions. However, the square and 
triangular shaped orifices were more effective mixing configurations 
at lower j conditions. 

Nomenclature 


Variables 


AMIX 

= 

overall enthalpy variance from h e q, Eq 4 

AHOT 

= 

enthalpy variance above hgq, Eq 5 

ACOLD 

= 

enthalpy variance below heq, Eq 6 

B 


area determined half-width of the 
distribution function. Figure 2 

C 


arithmetic average of the distribution 
function. Figure 2 

DF 

= 

distribution function (histogram). Figure 2 

DR 

= 

density ratio (jet/mainstream) 

DP/P 

= 

pressure drop across the mixing wall 

f 

= 

nondimensional specific enthalpy, Eq 2 

h 

= 

specific enthalpy 

J 

= 

momentum flux ratio (jet/mainstream) = M2/DR 

L 

= 

slot length (long dimension) 

rh 

= 

mass flowrate 

M 

= 

mass flux ratio (jet/mainstream) = DR Vj et /U ma j, 

MR 

= 

mass flowrate ratio (jet/ mainstream) 

R 

= 

radius of the mixing section 

S 

= 

spacing between orifice centers 

T 

= 

temperature — K 

Umain 

= 

inlet mainstream axial velocity 

Vjet 

= 

radial velocity of the jet 

VR 

= 

velocity ratio (jet/mainstream) = Vj e t/U ma in 

W 

= 

slot width (short dimension) 

X 

= 

axial distance from the leading edge of the 
orifice 

$ 


slot slant angle with respect to the axial 
direction 
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Introduction 

A number of advanced technologies are being developed for 
incorporation in the future generation turbopropulsion engines to 
simultaneously achieve challenging goals of improved efficiency, 
reduced life-cycle cost, lower noise, and exhaust gaseous and partic- 
ulate emissions. The major exhaust species of concern arc unburned 
hydrocarbons (UHC), carbon monoxide (CO), smoke, and oxides of 
nitrogen (NO x ). The UHC, CO, and smoke emissions from modem 
gas turbines have been reduced significantly by making design 
improvements (e.g., good fuel spray atomization quality and mini- 
mum quenching during low power operation) with little impact on 
life-cycle cost and safety. Although some progress has been made in 
lowering the NO x emissions by reducing the combustor length and 
the areas of stoichiometric combustion, the limitations of conven- 
tional single-stage turbopropulsion combustors are well known in 
regard to how much NOx emissions can be reduced without impact- 
ing operational safety. 

A promising approach to achieve significant NOx reductions for 
future-generation engines with very high cycle pressure ratios and 
turbine inlet temperatures is staged combustion (Rizk and Mongia, 
1990, 1991). The first stage is a rich bum zone in which fuel reacts at 
temperatures significantly below the stoichiometric temperature to 
prevent the formation of NO x . Typical rich burn zones operate at 
equivalence ratios between 1.3 to 1.8 and reaction temperatures of 
about 1800 K. The second zone is a quick mixing region in which the 
hot gases exiting the rich bum region are rapidly mixed with addi- 
tional air to minimize formation of NOx- The third zone is a lean 
burn zone typically operating at an equivalence ratio of about 0.7 to 
0.5. In this zone, the combustion process is completed by oxidizing 
the partially burned combustion products of the rich zone. 

As discussed by Rizk and Mongia (1990, 1991 ), RQL NOx emis- 
sions can be reduced by improving the mixing in the rich zone; 
increasing the rich zone length and residence time; using an opti- 
mum rich zone equivalence ratio; and achieving the quickest pos- 
sible rate of mixing in the quench zone. The latter has been shown to 
be the most critical process for controlling NOx emissions from RQL 
combustors (Rizk and Mongia, 1991). For example, NO x produced 
in the quench zone can be as high as 90% of the total NO x emissions 
of a current technology RQL combustor. 

To understand the fundamental processes within the mixing 
zone, a systematic detailed experimental and analytic effort is being 
conducted by several investigators (Hatch, Sowa, Samuelsen, 
Holdeman, 1992; Talpallikar, Smith, Lai, Holdeman, 1991; Vranos, 
Liscinsky, True, Holdeman, 1991; and Smith, Talpallikar, Holdeman, 
1991). Multidimensional gas turbine combustor models (e.g., COM- 
3D) can be used to analytically screen the effects of important design 
variables (e.g., AP, MR, rich zone equivalence ratio, J), geometric 
variables including shape, orientation, and number of air injection 
orifices on the quench zone mixing processes, and attendant NOx 
formation. This paper summarizes recent work in modeling 
nonreacting flow fields in a simple mixing zone test section for 
which experimental data have recently been reported by Hatch et al., 
1992. Parameters such as orifice geometry, orifice orientation, and 
number of orifices (equally spaced) around the circumferential 
direction were analyzed. Results are presented describing the 
relation of several parameters. This study also describes the 
interaction of multiple jets and how this affects the mixing process in 
a cylindrical duct. Acceptable agreement has been found between 
the qualitative and quantitative results using several statistical 


correlations to examine the numerical results to evaluate the mixing 
performance of a given configuration. 

Mathematical Model 

A production 3-D combustor code (COM-3D) solves the turbu- 
lent reading flow transport equations using the SIMPLE algorithm of 
Patankar and Spalding (Patankar, 1980). This program simulates 
turbulence by the two-equation k-e model (Launder and Spalding, 
1974) and combustion following vaporization is determined by a 
four-step chemical reaction model based on Arrhenius and eddy 
breakup concepts (Mongia and Reynolds, 1979). The transport 
equations for all dependent variables are of the following form as 
shown in Eq 1: 

div [pu4-(|l e ff /P r )grad(4)l = S^ (1) 

where p is the mixture density, u is the velocity, peff is the effective 
turbulent viscosity, P r is the effective Prandtl/Schmidt number, and 
S£ is the source term for the variable y. The following variables are 
computed by the three-dimensional code: (1) axial, radial, and swirl 
velocity components; (2) specific enthalpy and temperature; (3) tur- 
bulence kinetic energy and dissipation rate; (4) unburned fuel and 
composite fuel fraction; and (5) fuel spray trajectory and evaporation 
rate. 

The computational effort is significantly reduced by modeling a 
sector of the mixing section comprising a single orifice. Therefore, 
the shape of the sector was dependent on the number of orifices 
equally spaced in the circumferential direction. Periodic boundary 
conditions were applied in the circumferential direction. These 
conditions were necessary for the mixing configurations 


generating swirl during the mixing process generally characteristic 
of the slanted slot orifice configurations. No-slip and adiabatic 
boundary conditions were applied at the outer wall defining the 
inside wall of the quick mixing section tube. The computation at the 
center axis was extrapolated based on the nearest surrounding val- 
ues due to its singularity. No-gradient and slip boundary conditions 
were applied at the center axis. Fully developed profiles (no axial 
gradients) were assumed at the exit boundary condition. 

Geometric Configuration 

In this study, the quick mixing section was modeled as a constant 
diameter cylindrical "can" with a single row of equally spaced ori- 
fices. The outer wall diameter is 80 mm and the axial length of the 
mixing section extended from X/R = -1 to X/R = 5 where X is ref- 
erenced from the leading edge of the orifice. Sufficient axial distance 
was provided both upstream and downstream of the orifice to pre- 
vent any entry or exit effects from modifying the flow structure in 
the computational domain of interest which is X/R<2. The computa- 
tional grid of the domain typically consists of 20,000 to 30,000 control 
volumes. The grid was typically denser near the orifice to resolve 
the high velocity and temperature gradients resulting from the inlet 
of the crossflow jet. An orthogonal view of the computational grid is 
shown in Figure 1. The grid was also configured to allow smooth 
progressive volume change between adjacent control volumes to 
help speed up the convergence of the solution. 

The geometric configurations of the jet orifices are also shown in 
Figure 1 . Three different general shapes of orifices were modeled: 
square, slanted slots, and equilateral triangles, and a total of 21 dif- 
ferent configurations were analyzed. Multiple shapes, number, and 
orientations were analyzed and the different orifice configurations 
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Figure 1. Depiction of (a) geometric configuration of the mixing section; (b) orifice description; and 

(c) computational grid. 
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Table 1. 

Configuration description. 


CONF 

# 

ORIFICE 

TYPE 

SLANT 

ANGLE 

L/W 

# ORIFICES 
PER ROW 

BLOCKAGE 

ORIFICE 
T.E.* X/R 

1 

SQUARE 



6 

0-319 

0338 

2 

SOUARE 



8 

0382 

0309 

3 

SOUARE 



10 

0.405 

0.263 

4 

SOUARE 



12 

0.454 

0.245 

5 

SOUARE 



16 

0.520 

0.204 

6 

SLOT 

45 DEG 

1 

8 

0.525 

0.418 

7 

SLOT 

45 DEG 

2 

8 

0.432 

0383 

8 

SLOT 

0DEG 

4 

8 

0.189 

0.619 

9 

SLOT 

15 DEG 

4 

8 

0.300 

0.591 

10 

SLOT 

30 DEG 

4 

8 

0.521 

0.531 

11 

SLOT 

45 DEG 

4 

8 

0.565 

0.440 

12 

SLOT 

45 DEG 

6 

8 

0.684 

0.550 

13 

SLOT 

45 DEG 

8 

8 

0.750 

0.602 

14 

SLOT 

60 DEG 

4 

8 

0.671 

0.377 

15 

SLOT 

75 DEG 

4 

8 

0.735 

0.275 

16 

SLOT 

90 DEG 

4 

8 

0.765 

0.155 

17 

2-D SLOT 

90 DEG 

7 

8 

1.000 

0.106 

18 

TRIANGLE 

HASE DWNSTR 


8 

0.498 

0318 

19 

TRIANGLE 

BASE UPSTR 


6 

0.443 

0.451 

20 

TRIANGLE 

BASE UPSTR 


8 

0.478 

0387 

21 

TRIANGLE 

BASE UPSTR 


12 

0.633 

0.317 


NOTE: The effective area of the orifice divided by the area of the mixing section is 0.169 for all configura- 
tions in Table- 1. 

* T.E. is the trailing edge of the orifice 


are shown in Table 1. The effective orifice area (ACd) of the entire 
orifice row and orifice coefficient of discharge was maintained con- 
stant (849.5 nrim^ and 0.8, respectively) throughout the entire para- 
metric study. This also indicates that the area ratio between the 
orifice and the mixing section (orifice area/can area = 0.169) is also 
constant for all configurations described in Table 1. The model of the 
orifice perimeter was accomplished by a stair-stepping 
approximation since COM-3D is not a body-conforming code. For 
example, the orifice models of configurations 18 and 20 (Table 1) are 
slightly different due to computational grid density. Therefore, to 
increase the modeling accuracy, the orifice was typically defined 
with 30 to 100 control volumes. 

The mixing section was modeled at atmospheric pressure. The 
mainstream flow and jet flow conditions are shown in Table 2. The 
mainstream temperature, mainstream flowrate, and jet temperature 
were maintained constant throughout the parametric study. The 
variation in J was chosen as a representative range for typical gas 
turbine combustion systems. The mainstream flow at the inlet of the 
modeled mixing region was characterized by a uniform temperature 
and axial plug flow velocity profiles in the radial and circumferential 
plane. The air jet flow was characterized by a radial , uniform flow 
across the orifice effective area. The turbulence kinetic energy of the 
mainstream and jet flows were 0.3% of the mean velocities. The 
turbulent length scales of the mainstream flow were 2% of the can 
diameter, and the turbulent length scale of the jet was of the order or 
the orifice diameter. The inlet conditions for all orifices in the mixing 
section were equal to create a symmetrical input condition about the 
circumferential direction that was necessary for the sector analysis. 
The operating pressure and temperature range was chosen to 
simplify the experimental setup for future numerical validation. 


A typical numerical solution took about 250 iterations for full 
convergence with overall mass flow residuals of 0.05% of the total 
mixing section mass flowrate. All solutions were obtained using the 
Cray Y-MP and a typical converged solution took about 1/2 hr of 
CPU time. 

Results and Discussions 

The mixing performance for all configurations analyzed in this 
study were evaluated at X/R = 1. Two different methods were used 
to analyze the numerical results: 

1 . The results were analyzed qualitatively by visual observation of 
the enthalpy and velocity field solutions. The specific enthalpy 
plots presented in this paper describe a domain from X/R = -0.6 
to X/R = 2.4. In all axial-radial enthalpy plots, the air flows from 
left to right and the air flows into the paper for the radial-tangen- 
tial plots. The bottom boundary on all axial-radial plots consti- 
tutes the center axis of the mixing section and the upper bound- 
ary is the outer wall. The representative two-dimensional plane 
slices illustrated in this paper are: 

a. axial-radial plane slices through the center of the orifice (The 
axial direction is plotted in the abscissa.) 

b. radial-tangential plots at X/R = 1 downstream from the 
orifice 

The specific enthalpy results are presented as normalized values 
with respect to the overall heat differential between the main- 
stream flow and the jet flow. This normalized specific enthalpy 
is defined in Eq 2. 


Table 2. 

Operating conditions of the mixing section. 


CODE 

PREFIX 

T main 
degK 

V main 
m/sec 

T jet 
deg K 

V jet 
m/sec 

MR 

DR 

VR 

J 

DP/P 

% 

Equil 

enthalpy 

A 

588 

10.8 

300 

39.1 

0.96 

1.57 

3.61 

20.55 

0.708 

0.510 

B 

588 

10.8 

300 

52.1 

1.28 

1.57 

4.82 

36.55 

1.258 

0.439 

C 

588 

10.8 

300 

78.1 

1.92 

1.57 

7.23 

82.18 

2.831 

0343 


3 


h-h ipt 

f=r — dS- (2) 

n main n jet 

Note: This is a conserved scalar axially beyond the jet 
injection. 

The value of f varies from 0 to 1, where 0 is the value of the 
unmixed jet and 1 is the value of mainstream flow. Note that f = 
1 - 8, where 0 is as defined previously (Holdeman, 1991). The 
value of equilibrium specific enthalpy was calculated in an 
adiabatic system. The equilibrium enthalpy was also normalized 
in the same manner as the specific enthalpy. 

2. The numerical results were also quantitatively analyzed using 
two different statistical techniques. 

a. The mixing effectiveness was defined quantitatively by 
subdividing the specific enthalpy field at the radial-tangen- 
tial plane corresponding to X/R = 1 into 20 bins covering the 
entire enthalpy range from lowest enthalpy (jet) to the high- 
est (mainstream). The enthalpy of each computational grid 
control volume was calculated and sorted into its appro- 
priate bin, and this procedure was done with all control vol- 
umes. The definition of the distribution function (DF) is of 
the form as described in Eq 3 and Figure 2. 

*>L 

T, (mjjkhjjk) 

DF histogram = (3) 

hL+Ah 

With this information, a mass flux weighted specific 
enthalpy distribution density function (DF) was constructed 
from the data. The partial enthalpy corresponding to each 
enthalpy bin was then normalized with respect to the 
enthalpy contained at the plane located at X/R = 1. There- 


fore, ideal mixing would occur at a location downstream of 
the injection plane where the value of the corresponding dis- 
tribution function is 1.0 at the appropriate enthalpy bin, 
which coincides with the equilibrium enthalpy (heq). Sub- 
sequently, the mixing effectiveness was quantified by the 
shape of the histogram (DF) formed by plotting the normal- 
ized enthalpy bin values (abscissa) and distribution function 
DF (ordinate). To further quantify the shape of the distribu- 
tion function, the following parameters shown in Figure 2 
were defined: 

i. B(+) and B(-) are the "area determined" distribution half 
width. B(+) is the half width above the equilibrium 
enthalpy, which is the deviation from heq such that the 
area under the histogram (DF) from heq to B(+) is 1/2 
that of the overall distribution area above heq. The same 
is applicable for B(-) with the area below heq. The defi- 
nition of both B(+) and B(-) are shown pictoriaily in 
Figure 2. 

ii. C(+) and C(-) are the distribution half width area aver- 
age values. C(+) is the overall arithmetic average of the 
distribution function for the values above heq. The same 
is applicable to C(-). The definition of both C(+) and C(-) 
are shown pictoriaily in Figure 2. The area described by 
the rectangle with base of 2B(+) and height of C(+) 
equals the area of the DF histogram above heq. 

Similarly, the rectangle described by base of 2B(-) and 
height of C(-) equals the overall area of the histogram 
below heq. 

b. The performance of the mixing section was also evaluated by 
using statistical area weighted variance computations. The 
area weighted RMS values were calculated at X/R = 1. 

Three different parameters (mixing deviations) — AMIX, 
AHOT, and ACOLD — are described in Eqs 4 through 6. 


AMIX = 


i 



(4) 


0 



Ah 2 Ii (-) B (-) h C q B(+)2B(+) 

Figure 2. Graphical definition of the distribution function (DF), B(+), B(-), C(+), and C(-). 
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AHOT = 


(5) 


where 


1 


" 


21 

1 Yai 

1 jk * eq 


A 2^ A ) k 
A, °‘ jk 

T'main — ^jet 



f° r Tjk 


>T, 


eq 


C = 2.50 


ACOLD = 


■ 

_ 


1 Va. 

Tjk ~ Tgq 


i 

>1 

o 

*5r 

Tmain “ Tj e t 



f °r T )k < T txJ (6) 


From the definitions described in Figure 2, the best mixing is 
attained with the highest possible C(+) and C(-) and lowest 
possible B(+) and B(-). Similarly, the smallest deviation 
AMIX, AHOT, and ACOLD indicate the best mixing 
configuration. Note that AMIX is the square root of the sum 
of the squares of AHOT and ACOLD. The results described 
below are presented separately for each condition listed in 
Table 2. In this way, all geometric configurations with the 
same J also had the same MR and DR. The following four 
different parametric studies were performed: 

• Effect of number of orifices. Quick mixing configurations 
with square orifices were analyzed with a single row of 6, 
8, 10, 12, and 16 orifices. In addition, orifice configura- 
tions with the equilateral triangle (base upstream) were 
analyzed with 6, 8, and 12 orifices/row only at J = 20.55. 
Higher J conditions were not analyzed in this study. 

• Effect of slot aspect ratio. Five different slot aspect ratios 
were analyzed with a slot slant angle of 45 deg. The dif- 
ferent slot aspect ratios are: 1, 2, 4, 6, and 8. All mixing 
region configurations had a single row of 8 orifices. 

• Effect of slot slant angle. Seven different slanted slot ori- 
entations were analyzed with a slot aspect ratio of 4. The 
different slot slant angles are 0, 15, 30, 45, 60, 75, and 90 
deg with respect to the center axis of the quick mixing 
region. All mixing configurations had a single row of 8 
orifices. 

• Effect of orifice shape. Seven different orifice shape con- 
figurations were analyzed: square, equilateral triangle 
(base downstream), equilateral triangle (base upstream), 
diamond, flow aligned slot, transverse slot, and continu- 
ous (2-D) slot. All mixing region configurations had a 
single row of orifices, and the overall row ACd for all 
configurations was maintained constant at 849.5 mm^. 

Results for I - 20.55 

Effect of Number of Orifices on Mixing 

The results shown in Figure 3 indicate a significant decrease in 
jet penetration with increase of the number of orifices. The 
configurations with 12 orifices/row had a larger blockage fraction 
(0.402) compared to 6 orifices/row (0.284). The blockage fraction of 
the jet orifice (Table 1) is defined as the circumferential projection of 
the orifice divided by the spacing between orifice centers. The 
blockage fraction indicates a value of the mainstream flow blockage 
produced by the jet at the location where the jet enters the mixing 
section. A higher blockage would effectively prevent the hot 
mainstream flow from occupying the volume immediately 
downstream of the orifice. Since the structure of the mixed flow in a 
cylindrical mixing section is highly area dependent, it is generally 
advantageous to have an orifice configuration with a higher blockage 
fraction to produce a more uniformly mixed flow structure toward 
the outer section of the mixing section. The enthalpy plot solutions 
indicate the coolest region behind the orifice for the configuration 
with 12 orifices/row, therefore, being the best mixing configuration. 
The statistical parameters also convey the same conclusion as the 
qualitative analysis as shown in Figure 8a. According to Holdeman 
(1991), the optimum number of orifices (n) in a can is shown in Eq 7. 
The optimum number of orifices calculated using Eq 7 for ] = 20.55 
is 8. 


n = 


*V2l 

C 


(7) 


The jet penetration produced with the configuration with 8 ori- 
fices/row is significantly different compared to that of 6 orificcs/row 
due to the change in jet flow structure near the center core of the 
mixing section. As the number of orifices increases, individual jets 
merge into a single structure that interacts differently with the main- 
stream flow. Although the jet penetration from a large number of 
smaller orifices is less than a smaller number of larger orifices, typi- 
cally, configurations producing a merged jet structure have consid- 
erably less penetration compared to the configurations producing 
more distinct isolated jets. The jet penetration is affected by the high 
static pressure at the center core, especially when the hot core is vir- 
tually shielded from the rest of the mixing section due to the merged 
jet structure. This therefore inhibits the penetration of the merged 
jets as shown in Figure 3 for configurations with 8, 10, and 12 ori- 
fices/row. The statistical parameters indicate a significant increase 
in the mixing performance with the increase in the number of ori- 
fices/row, especially between 6 and 8 orifices/row. Both area and 
mass flux weighted parameters indicate similar trends. 

The results for the equilateral triangles in Figure 4 seem to 
indicate a slight decrease in mixing performance with increasing 
number of orifices due to the rapid decrease in jet penetration. This 
trend, however, is opposite to the trend obtained with the square 
orifices. Based on these analytical results, the optimized triangular 
configuration at J = 20.55 has 8 orifices/row. 

Effect of Slot Aspect Ratio 

The results shown in Figure 5 and 8b indicate little improvement 
in the mixing performance with increase in slot aspect ratio from 
L/W from 1 to 6. Also note that the axial distance of injection 
increases with aspect ratio as shown in Table 1. A significant 
improvement was obtained by increasing the L/W to 8, however 
both area and mass flux weighted parameters indicate similar 
results. 

Effect of Slot Slant Angle on Mixing 

The results shown in Figures 6 and 8c indicate a slight 
improvement in the mixing performance with increase in slot slant 
angle. The mass flux weighted parameters indicate a higher 
sensitivity to the change in the slant angle compared to the area 
weighted parameters. The mass flux weighted parameters are 
significantly influenced by the value of the axial velocity distribution 
at plane X/R = 1 and these results are influenced by the amount of 
swirl produced in the mixed flow structure. Furthermore, the 
amount of swirl imparted to the mixed flow depends on the slant 
angle. Highest swirl was produced with the 45-deg slant angle slot. 
Both 0-deg and 90-dcg slanted slot configurations produced no swirl 
(as evidenced by the symmetric distribution in Figure 6). The 60-deg 
slant angle slot produced the best mixing at J = 20.55. 

Effect of Orifice Shape 

The results shown in Figure 7 indicate that optimum mixing per- 
formance is obtained with the configurations with the moderate jet 
penetration which are produced with both the square and triangle 
(base upstream) orifices. Excessively high jet penetration (such as 
the flow aligned slot) or low (such as the continuous [2-D] slot) pro- 
duce nonoptimal cases. The statistical results are shown in Table 3 

Table 3. 

Effect of orifice shape on mixing for J = 20.55, MR = 0.96, and DR = 1.57. 

TAIII.K-.V F.I- I-KCT OK ORIKICK SHARK ON MIXING KOR J 211.55, MR-0 96. AND UR-1. 5" 


ORIFICE CONTOUR 

!*{♦) 

CO) 

B(-) 

C(-) 

AMIX 

AHOT 

ACOLD 

SOI ARE 

0.119 

0.101 

0.098 

0.151 

0.157 

0.14 

0.07 

TRN« RASE UPST 

0.237 

0.06 

0.165 

0.067 

0.213 

0.188 

0.099 

TRNCI BASF DWNST 

0.124 

0.088 

0.089 

0.16 

0.129 

0.15 

0.076 

DIAMOND 

0.18 

0.073 

0.165 

0.072 

0.19 

0.156 

0.108 

SLOT 0=0 deg 

0.215 

0.063 

0.165 

0.07 

0.202 

0.164 

0.1 IX 

SLOT 0=90 doc 

0.155 

0.066 

0.09 

0.162 

0.17 

0.135 

0.103 

:d si.ot 

0.395 

0.021 

0.325 

0.031 

0.391 

0JI3 

0.235 
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10 ORIFICES / ROW 12 ORIFICES / ROW 

a. Axial-radial plane 
through orifice center 





b. Radial-tangential plane 
at X/R=1 


f 



1.00 

0.93 

0.86 

0.79 

0.71 

0.64 

0.57 

0.50 

0.43 

0.36 
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Figure 3. Effect of the number of square orifices on mixing for J = 20.55, MR = 0.96, DR = 1.57, 

hgq = 0.510. 
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Figure 4. Effect of the number of triangular orifices on mixing for J = 20.55, MR = 0.96, DR = 1 .57, 

heq = 0.510. 
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L/W = 8 


a. Axial-radial plane 
through orifice center 



b. Radial-tangential plane 
at X/R=1 


Figure 5. Effect of slot aspect ratio (L/W) on mixing for J = 20.55, MR = 0.96, DR = 1 57, slang angle = 

45 deg, hpq = 0.510. 
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b. Radial-tangential plane 
at X/R=1 


Figure 6. Effect of slot slant angle on mixing for J = 20.55, MR = 0.96, DR = 1.57, L/W = 4, heq = 0.510. 
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Figure 7. Effect of orifice shape on mixing for J = 20.55, MR = 0.96, DR = 1.57, hgq = 0.510. 
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Figure 8. Statistical results of the parametric study at J = 20.55, MR = 0.96, and DR = 1 .57. 
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Results for I = 36.50 

Effect of Number of Orifices 

Both visual and statistical analyses shown in Figures 9 and 12a 
indicate an improvement in mixing with increase in the number of 
orifices. The increase in number of orifices increases the blockage 
fraction and reduces the mainstream hot gas recirculation behind the 
orifice. This produces a more uniform temperature distribution 
behind the mixing orifice at the expense of a hot unmixed center 
core. Both area weighted and mass flux weighted parameters indi- 
cate the same conclusion. The optimum number of orifices calcu- 
lated with Eq 7 yields 11 orifices/ row, which also agrees with the 
present results. 

Effect of Slot Aspect Ratio 

The results are shown in Figures 10 and 12b. These results show 
a significant improvement in mixing with increase in slot aspect 
ratio. The statistical results indicate an almost direct relationship 
between mixing improvement and increase in slot aspect ratio 
despite the fact that large L/W injects the mixing airflow over a 
longer axial distance. Optimum performance was obtained with the 
L/W = 8 slot. Both area weighted and mass flux weighted param- 
eters indicated the same conclusion. Since the statistical parameters 
are all highly area weighted, increased performance can be obtained 
by reducing the over-penetration of the L/W = 1 slot at the expense 
of producing a hot center core as shown for L/W = 8. However, the 
hot core in the center occupies a small portion (less than 5%) of the 
total area in the plane at X/R = 1, and this does not significantly 
affect the statistical results. 

Effect of Slot Slant Angle 

The results shown in Figures 1 1 and 12c indicate better mixing is 
produced with the larger slant angle configurations. As jet penetra- 
tion decreases with increase in slot slant angle, the area toward the 
outer wall of the mixing section becomes cooler and therefore 
increases the temperature uniformity. In addition, increase in slot 
slant angle up to 45 deg produces increased counterswirling flows, 
which stratify the temperature distribution in the radial direction 
and inhibit the entrainment of cold jet flow in the region behind the 
orifice. Some deviation was obtained between the area weighted and 
the mass flux weighted results. The mass flux weighted parameters 
are sensitive in areas where the axial velocity distribution is 
nonuniform, therefore, being highly dependent on the jet penetration 
region (where the velocity and gradients are high). The mass flux 
weighted parameters, unlike the simple area weighted averages, 
more accurately reflect the mixing tendency in areas where the flow 
gradients are high. Best mixing tendencies were obtained with either 
the 75- or 90-deg slant angle slots due to the best combination of 
optimum jet penetration, and low recirculation strength, and high 
blockage fraction. 

Results for I = 82.18 

Effect of Number of Orifices 

The results are shown in Figures 13 and 16a. The best perfor- 
mance was obtained with 16 orifices/ row at this high J mostly due to 
a significant reduction in the jet penetration. Configurations with 8- 
12 orifices overpenetrated due to the high J condition. These results 
are not surprising since results using Eq 7 indicate an optimum 
number of orifices of 16 at J = 82.18. 

Effect of Slot Aspect Ratio 

The results are shown in Figures 14 and 16b. The results indicate 
an optimum mixing performance with the large aspect ratio slot 
L/W = 8. The increase in slot aspect ratio decreases the jet penetra- 
tion as shown previously for the lower J conditions. Overpenetration 
is produced with the slot with aspect ratio of L/W = 4 and a high 
temperature nonuniformity is also observed due to the high temper- 
ature entrained behind the orifice. The area weighted parameters. 


which are nonvelocity dependent, show a different trend compared 
to the mass weighted parameters. This is due to the high velocity 
gradients produced by the swirling flows which are highest for the 
aspect ratio orifice (L/W = 4). After incorporating the velocity effect 
into the area weighted parameters, the results were similar in trend 
to the distribution function mass weighted parameters. 

Effect of the Slot Slant Angle 

The results are shown in Figures 15 and 16c. The results indicate 
a similar trend as described in the lower J configurations. Increased 
mixing performance was obtained with the large slant angle; how- 
ever, the optimum configuration producing best mixing was the con- 
figuration with the 75-deg slanted slot. Slanted slot configuration 
with up to 45-deg slant produced overpenetration at this high J; 
therefore, hot mainstream flow was entrained behind the orifice and 
produced an undesirable mixing structure. 

Conclusions 

1. The best orifice configuration depends significantly on the pene- 
tration depth of the jet to prevent the hot mainstream flow from 
being entrained behind the orifice. Orifices producing strong 
swirl in the mixing region have an adverse effect due to the tem- 
perature stratification toward the outer wall. 

2. The simple area weighted statistical parameters correlate well 
with the mass flux weighted parameters in almost all analyzed 
configurations. Differences were found however in flow struc- 
tures where swirls were high and axial velocity distributions 
contain high gradients. In such cases, the area weighted param- 
eters are somewhat deficient in quantifying the mixing tendency. 
The flow structure in a circular mixing section is highly weighted 
toward the outer wall and any flow structure affecting this area 
significantly affects the overall results. 

3. The increase in the number of orifices per row increases the mix- 
ing performance for low J due to the increase in the orifice block- 
age fraction. Significant mixing improvement was also 
attributed to the merging tendency of the individual jets in a cir- 
cular domain. Jet overpenetration produced nonoptimal mixing 
configurations. 

4. Higher slot slant angles and aspect ratios are generally the best 
mixing configurations at higher J conditions. However, at the 
lowest J conditions analyzed in this study, the square and trian- 
gular orifices have better mixing. Therefore the design of the 
quick mixing section strongly depends on the operating condi- 
tion of the combustion liner (DP/P and or J) and these conditions 
must be known a priori before selecting the appropriate orifice 
configuration for best mixing. 
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through orifice center 



b. Radial-tangential plane 
at X/R=1 

Figure 10. Effect of slot aspect ratio (L/W) on mixing for ] = 36.55, MR = 1.28, DR = 1.57, slant angle = 

45 deg, h e q = 0.439. 
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b. Radial-tangential plane 
at X/R=1 



Figure 11. Effect of slot slant angle on mixing for J = 36.55, MR = 1.28, DR = 1.57, L/W = 4, heq = 0.439. 
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Figure 12. Statistical results of the parametric study at J = 36.55, MR = 1.28, and DR = 1.57. 
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a. Axial-radial plane 
through orifice center 


b. Radial-tangential plane 
at X/R=1 




Figure 13. Effect of the number of square orifices on mixing for J = 82.18, MR = 1.92, DR = 1.57, 

heq = 0.343. 
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Figure 14. Effect of slot aspect ratio (L/W) on mixing for J = 82.18, MR = 1 .92, DR = 1 57 , slant angle = 

45 deg, heq = 0.343. 
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b. Radial-tangential plane 
at X/R=1 


Figure 15. Effect of slot slant angle on mixing for J = 82.18, MR = 1.92, DR = 1.57, L/W = 4, heq - 0.343. 
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c. 
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Figure 16. Statistical results of the parametric study at f = 82.18, MR = 1.92 and 1 57 
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